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MASS AND DUST IN THE DISK OF A SPIRAL LENS GALAXY 
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ABSTRACT 

Gravitational lensing is a potentially important probe of spiral galaxy structure, but only a few cases 
of lensing by spiral galaxies are known. We present Hubble Space Telescope and Magellan observations 
of the two-image quasar PMN J2004-1349, revealing that the lens galaxy is a spiral galaxy. One of the 
quasar images passes through a spiral arm of the galaxy and suffers 3 magnitudes of y-band extinction. 
Using simple lens models, we show that the mass quadrupole is well-aligned with the observed galaxy 
disk. A more detailed model with components representing the bulge and disk gives a bulge-to-disk mass 
ratio of 0.16 ± 0.05. The addition of a spherical dark halo, tailored to produce an overall flat rotation 
curve, does not change this conclusion. 

Subject headings: galaxies: spiral, structure, halos — dark matter — gravitational lensing — dust, 
extinction — quasars, individual (PMN J2004-1349) 



1. INTRODUCTION 

A general property of spiral galaxies is that they have 
three distinct mass components: a bulge, a disk, and an 
extended halo of dark matter. The bulge and disk can be 
detected and measured in optical images. The evidence for 
dark matter comes from observations of luminous tracers, 
such as the dynamics of galactic satellites, tidal tails, and 
globular clusters; and, most famously, H I and stellar ro- 
tation curves (as reviewed recently by Sofue & Rubin 2001 
and Combes 2002) . Instead of exhibiting a Keplerian fall- 
off, the rotation curves are flat, or nearly so, at a distance 
of many optical scale lengths from the galaxy center. 

Where the galaxy light is still appreciable, it is not clear 
what fraction of the total gravitational force that produces 
the rotation curve is provided by the luminous matter, or, 
equivalently, what is the shape and density profile of the 
dark halo. Even the extreme hypothesis that the inner 
rotation curve is produced by the luminous disk alone (a 
"maximum disk" ; Van Albada & Sancisi 1986) is still the 
subject of debate (see, e.g., Sackett 1997, Corteau & Rix 
1999, Palunas & Williams 2000). Dark haloes are tradi- 
tionally visualized as spherical, although there is evidence 
that they are flattened in the direction perpendicular to 
the disk, with density distributions that have a three- 
dimensional axis ratio ~0.5 on scales of ^^15 kiloparsecs 
(see Sackett 1999 and references therein). 

Gravitational lensing has long been recognized as a valu- 
able tool for studying the mass distribution of galaxies, 
and of spiral galaxies in particular (see, e.g., Mailer, Flo- 
res, & Primack 1997, Keeton & Kochanek 1998, Koop- 
mans, de Bruyn, & Jackson 1998, MoUer & Blain 1998). 
The image positions and magnifications of a strongly 
lensed quasar depend on the projected mass within the 
radius bounded by the quasar images, which is typically 
comparable to the optical radius of the lens galaxy. Fur- 
thermore, gravitational deflection depends on the total in- 



tervening mass, regardless of its luminosity or internal dy- 
namics. 

Unfortunately, few spiral galaxy lenses are known. Of 
the present sample of approximately 80 multiple-image 
quasars, only 4 are confidently known to be produced by 
spiral lenses. Selection effects favor the discovery of el- 
liptical lenses over spiral lenses: spirals have a smaller 
multiple-image cross section, spirals produce multiple im- 
ages with smaller angular separations, and spirals contain 
large amounts of dust that can extinguish one or more 
images (see, e.g., Wang & Turner 1997; Perna, Loeb, & 
Bartelmann 1997; Keeton & Kochanek 1998; Bartelmann 
& Loeb 1998). 

Two of the 4 spiral lenses described previously are 
poorly suited for mass modeling. The spiral lens galaxy 
of B0218-F357 (Patnaik et al. 1993) is so crowded by the 
bright quasar images that even its position has not yet 
been measured accurately, despite imaging with the Hub- 
ble Space Telescope (Lehar et al. 2000). Likewise, there 
appears to be a spiral lens in PKS 1830-211 (Pramesh Rao 
& Subrahmanyan 1988) but the field is so crowded with 
stars that two very different conclusions about the lensing 
scenario have been drawn from the same data (Courbin et 
al. 2002, Winn et al. 2002). 

The "Einstein cross" Q2237-I-0305 is produced by a spi- 
ral galaxy (Huchra et al. 1985), and is also unusual in an- 
other way: the galaxy has the lowest redshift (z = 0.0394) 
of all known lens galaxies. This causes the quasar images 
to appear very close to the galaxy center (<1 kpc), where 
they are sensitive mainly to the mass in the bulge rather 
than the disk or halo. Nevertheless, some interesting in- 
formation about mass on larger scales has been obtained. 
For example, Schmidt, Webster, & Lewis (1998) estimated 
the mass of the bar that is seen in optical images, from its 
shearing effect on the image configuration. Trott & Web- 
ster (2002) argued that the disk is sub-maximum using 
model constraints from both the quasar image configura- 



^ Harvard-Smithsonian Center for Astrophysics, 60 Garden St., Cambridge, MA 02138 
^ National Science Foundation Astronomy & Astrophysics Postdoctoral Fellow 
^ Princeton University Observatory, Princeton, NJ 08544 

* Departamento de Astronomfa y Astrofi'sica, Facultad de Ffsica, Pontificia Universidad Catolica de Chile, Casilla 306, Santiago 22, Chile 
^ Department of Physics, Massachusetts Institute of Technology, 77 Mass. ave., Cambridge, MA 02139 

1 



2 



Winn, Hall, & Schechter 



tion and H I rotation measurements at larger galactocen- 
tric distances. 

The other well-studied case of spiral lensing is 
B1600+434 (Jackson et al. 1995, Jaunsen fc Hjorth 1997), 
in which two quasar images bracket the nearly edge-on disk 
of the lens galaxy. Mailer, Flores, & Primack (1997) pre- 
sented the first models for this system consisting of both a 
halo and a disk, which were extended by Mailer et al. 2000 
to determine the allowed combinations of disk mass and 
halo ellipticity. Koopmans, de Bruyn, & Jackson (1998) 
achieved a similar goal using models with disk, bulge, and 
halo components. However, this system has two undesir- 
able properties for mass modeling: the lens galaxy has a 
central dust lane that makes its position hard to measure 
accurately, and there is a massive neighboring galaxy that 
adds complexity and uncertainty to the models. In addi- 
tion, the quasar images are nearly collinear with the galaxy 
center, an accidental symmetry that makes it difficult to 
test whether the overall mass distribution is aligned with 
the galaxy disk; previous studies have assumed this is the 
case. 

In this paper we describe a spiral lens system that of- 
fers a new opportunity for mass modeling. The system, 
PMN J2004-1349, is a two-image quasar originally dis- 
covered in a radio lens survey (Winn et al. 2001). The 
lens galaxy was identified in ground-based optical images 
but the angular resolution of those images was not good 
enough to determine the galaxy's position or morphology. 
The quasar and lens galaxy redshifts are unknown. In 
§ 2 we present Hubble Space Telescope (HST) optical im- 
ages showing that the lens galaxy is a spiral galaxy. The 
quasar images have very different optical colors. In § 3 
we present Magellan data that extend the color measure- 
ments to near-infrared wavelengths, and we test whether 
the color differences are consistent with differential extinc- 
tion by dust in the lens galaxy. 

In § 4 we present mass models. The lens galaxy does not 
have obviously massive neighbors, its position is known 
with greater accuracy than the problematic cases men- 
tioned above, and it does not lie along the line between 
the quasar images. Yet despite these advantages, it is still 
only a two-image system, and two-image systems do not 
provide many constraints on lens models. Our approach is 
to determine what can be learned from the simplest plau- 
sible mass models, and then consider more complicated 
models. In § 4.1 we use only the quasar image positions 
and fluxes to test whether the mass quadrupole is aligned 
with the luminous disk, since this is the first system for 
which such a test has been possible. Then in § 4.2, we use 
the HST surface photometry to constrain a bulge+disk 
model, in order to measure the bulge-to-disk mass ratio. 
In both §§4.1 and 4.2, we briefly consider the implications 
of an additional assumption: a flat rotation curve. Finally, 
in § 5 we summarize and discuss future observations that 
would provide more constraints on spiral galaxy structure. 



Parameter 




Value 


R.A.NE- R.A. 


■ SW 


981.47 ± 1 mas 


decl.NE— decl. 


SW 


552.34 ± 1 mas 


R.A.NE- R.A. 


■Gal 


780 ± 5 mas 


decl.NE— docl. 




247 ± 5 mas 


^NE/'/Jri\V 




1.0 ± O.i 



Table 1 — Lens model constraints 



2. HST OBSERVATIONS 

On 2001 May 21, we observed J2004-1349 with the 
HST^ and WFPC2 camera (program id 9133). Subse- 
quently, on 2001 June 24, the same fleld was observed by 
Beckwith et al. (id 9267), for a high-redshift supernova 
search. All together there were 20 dithered exposures us- 
ing the F814W filter (hereafter, /), totaling 17.2 ksec; and 
4 dithered exposures using the F555W filter (V), total- 
ing 5.1 ksec. The target was centered in the PC chip, 
which has a pixel scale of 0'.'0456 pixel^^. The exposures 
were combined and cosmic rays rejected using the driz- 
zle package (Fruchter & Hook 2001) and other standard 
IRAF^ procedures. 

Figure 1 shows the final /-band image, at low contrast 
(left panel) and high contrast (right panel). The quasars 
are labeled NE and SW. The lens galaxy, G, is a spiral 
galaxy. It has a prominent bulge and a disk extending 
nearly north-south. Two loosely wound spiral arms are 
visible, spiralling clockwise as they emerge from the east- 
ern and western sides of the bulge and extend north and 
south. The spiral arms are more obvious in the image 
shown in Figure 3, in which the quasars and an elliptically 
symmetric galaxy model have been subtracted using a pro- 
cedure described below. The SW quasar passes through 
the southern spiral arm, close to the bulge. The /-band 
images also revealed a faint neighboring galaxy, G2, and 
another object, X. Object X seems extended but it is too 
faint to rule out the possibility that it is a foreground star. 
In V-band, the quasars and the lens galaxy bulge were de- 
tected, but with a lower signal-to-noise ratio. From the 
HST images we obtained the following information: 

Relative positions. We measured the relative positions 
of the quasars and lens galaxy with point-spread-function 
(PSF) fitting, using a nearby star as an empirical PSF 
(star #3 from the finding chart of Winn et al. 2001). We 
modeled the quasars as point sources, and the galaxy bulge 
as an elliptical exponential profile convolved with the PSF. 
Errors were estimated by the variance in results of fits to 
three sub-images, each of which was constructed from only 
one-third of the data. The quasar separation agreed with 
the more precise value determined by Winn et al. (2001) 
using vc;ry-long-baseline interferometry. The lens galaxy 
position is given in Table 1 along with other constraints 
used in lens models (§ 4). The result did not change signifi- 
cantly when the lens galaxy was modeled with an elliptical 
Gaussian profile or a de Vaucouleurs profile. 

Flux ratio. The quasar flux ratio (NE/SW) is 8:1 in 
/-band and 20:1 in F-band, as compared to 1:1 at radio 
wavelengths. The SW quasar image is signiflcantly redder 
than the NE quasar image. In § 3 we describe data that ex- 

^ Data from the nasa/esa Hubble Space Telescope (HST) were obtained at the Space Telescope Science Institute, which is operated by the 
Association of Universities for Research in Astronomy, Inc., under NASA contract NAS 5-26555. 

The Image Reduction and Analysis Facility (iraf) is a software package developed and distributed by the National Optical Astronomical 
Observatory, which is operated by aura, under cooperative agreement with the National Science Foundation. 
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Fig. 1. — HST Image of J2004-1349. North is up and east is left. The quasars are separated by 1''13. 



tend the measurements to near-infrared wavelengths and 
discuss the interpretation. 

Position angle and axis ratio of the disk. By overplot- 

ting ellipses on the /-band image, we visually measured 
the position angle of the galaxy disk, 9i = — 6?5 ± 3?5, 
and its projected axis ratio, qi — 0.41 ± 0.05. We use the 
subscript I to distinguish properties of the luminosity dis- 
tribution from properties of the mass distribution, which 
we will denote with the subscript m. Assuming the galaxy 
is intrinsically a flat circular disk, the inclination is given 
hyi = cos^^qi = 66° ± 3°. 

Standard magnitudes. To measure the galaxy flux, we 
used a I'.'S x 2f.'5 rectangular aperture centered on the 
bulge, after subtracting the quasars. Table 2 gives the 
magnitudes, using the Dolphin (2000) zero points^ of 
21.654 for I and 22.551 for V. The total magnitudes agree 
with the ground-based photometry of Winn et al. (2001), 
but the magnitudes of the individual components do not 
agree. This is because of the much poorer angular res- 
olution of the ground-based images. In particular, what 
was previously identified as the SW quasar in the ground- 
based images is now known to be mainly light from the 
lens galaxy bulge. 



Component 


F814W 


^ / 


F555W 


^ V 


NE 


22.06 ± 


0.02 


24.54 ± 


0.04 


SW 


24.56 ± 


0.11 


27.79 ± 


0.41 


G 


21.62 ± 


0.03 


24.30 ± 


0.16 



Table 2 — HST photometry. Error estimates represent sta- 
tistical error only, and do not include the zero point error of 
wO.05 mag or GTE effects. 



Surface brightness profile. Figure 2 shows the surface 
brightness of the galaxy averaged over elliptical contours. 
We experimented with several parametric fits to the sur- 
face brightness profile, and found the best results for a 
model consisting of the sum of two exponential profiles 
with scale lengths 0'.'064 and 0'.'65.^ The solid line shows 
the profile of the model after convolving with the PSF. 
The dashed and dotted lines show the contributions from 
each of the two components. The model is a good fit ex- 
cept for the "bump" at a semi-major axis of O'.'S due to 
the spiral arms. The arms are obvious in Figure 3, which 
shows the residual image after subtracting the quasars and 
galaxy model. In this model, the bulgc-to-disk flux ratio is 
0.33 ± 0.07, where the disk flux includes the total residual 
fiux within the rectangular aperture described above. 



3. DIFFERENTIAL EXTINCTION 

The HST images showed that the SW quasar image is 
redder than the NE quasar image. This trend was con- 
firmed with near-infrared images obtained on 2002 June 29 
with the Baade 6.5m (Magellan I) telescope and Classic- 
Cam, a 256^ HgCdTc imager with 0'.'115 pixels. We ob- 
tained 52.5 minutes' integration in Jg and 44 in Kg, both 
in 0'.'4 seeing and non-photometric conditions. The data 
were reduced in the standard fashion. We measured the 
quasar flux ratio using the same PSF-fitting procedure 
that was used for the HST images, except that in this case 
we locked the relative positions (and structural parame- 
ters of the galaxy) at the HST-derived values and allowed 
only the fluxes to vary. The uncertainty was estimated by 
the spread in results obtained for different choices of the 



^ As updated at http://www.noao.edu/staff/dolphin/wfpc2_calib/ and corrected to infinite aperture. 

^ We note that this decomposition is not unique; for example, a poorer but still reasonable fit can also be achieved with a de Vaucouleurs bulge 
(iieff = 0^'2) instead of an exponential bulge. 
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Fig. 2. — (left) Surface brightness of the lens galaxy, averaged over ellipses with position angle 6i = — 6?5 and axis ratio qi = 0.41. Error bars 
show the variance of the points used to compute the average. The lines show the corresponding surface brightness profile of the components 
in the model image. 



Fig. 3. — (right) Residual /-band image (= data — quasar model — galaxy model). 



lens galaxy profile. Table 3 gives the quasar flux ratio at 
optical and near-infrared wavelengths. 



Band 


Central wavelength (/im) 


Flux ratio (NE/SW) 


V 


0.555 


20 ±8 


I 


0.814 


8.3 ±0.7 


Js 


1.24 


2.7 ±0.4 


Ks 


2.16 


1.44 ±0.09 



Table 3 — Quasar flux ratio. 



Because gravitational lensing does not alter the wave- 
length of photons, the different colors of NE and SW re- 
quire explanation. The most likely explanation is that SW 
is being reddened by dust as its light passes through the 
spiral arm of the lens galaxy. Differential reddening has 
been observed in many lens systems (Nadeau et al. 1991, 
Falco et al. 1999), including the spiral lens B1600±434 
(Jaunsen & Hjorth 2002). The main alternative explana- 
tions are microlensing and intrinsic variability. Microlens- 
ing causes color changes if the angular size of the quasar 
is unresolved (which is always the case at optical wave- 
lengths) and varies with wavelength, because microlens- 
ing magnification depends on source size (Wambsganss & 
Paczynski 1991; see Wucknitz et al. 2003 for a recent ex- 
ample). Intrinsic variability can also cause color changes 
when the time delay is longer than the time scale of varia- 
tions and the degree of variability depends on wavelength. 
But neither of these alternative hypotheses is expected to 
produce color differences as large as observed here, and 
neither one would naturally explain why SW is systemat- 
ically redder than NE. 

To verify that the dust explanation is reasonable, we 



compare the wavelength-dependence of the quasar flux ra- 
tio, /a(NE)//a(SW), with what would be expected from 
dust reddening. For simplicity we assume that NE is unaf- 
fected by dust, and that intrinsic variability is negligible. 
Because the radio flux density ratio is 1:1, any differences 
in the observed magnitudes m\ should be due to extinction 

mA(SW) - mA(NE) = -2.51ogio = Aa(SW). 

We compare; the observations with an extinction law 
Ax/Av determined by Cardelh, Clayton, & Mathis (1989), 
using the standard Galactic value Ry = 3.1, and adjusting 
the wavelength scale as appropriate for the lens redshift. 
There are 2 adjustable parameters, E{B — V) and ziens) 
and 4 flux ratio measurements. The best fit (x^ = 1.4) is 
achieved for E{B — V) = 0.97 and ^lens = 0.16, which are 
both reasonable values. The measurements and the fitted 
extinction curve arc shown in Figure 4. We conclude that 
differential extinction is a good explanation for the color 
difference, although it is certainly possible that microlens- 
ing and intrinsic variability also contribute to the color 
difference. 

Following ,Jcan & Surdcj (1998), this can be regarded 
as a measurement of the lens redshift, but it is a crude 
one: the range of lens redshifts giving Ax^ < 1 is 0.03 < 
^lons < 0.36. In addition, the result depends on the choice 
of Rv = 3.1. Larger values of Rv favor smaller lens red- 
shifts and smaller values of Rv allow for large lens red- 
shifts. As an example, a solution with Ry = 2.1 and 
-^lens = 0.9 is also plotted in Figure 4. 
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Fig. 4. Alcasurcd extinction of quasar image SW (see Eq. 1) 
compared to two model extinction laws from Cardelli, Clayton, & 
Mathis (1989). 

4. LENS MODELS 

4.1. Power-law models 

We consider lens models constrained by the image po- 
sitions and magnification ratio (as approximated by the 
radio flux density ratio) given in Tabic 1. The uncertain- 
ties in the quasar positions and magnification ratio were 
enlarged from the observational uncertainties to 1 milliarc- 
second and 10%, respectively, to account for systematic 
effects due to mass substructure (see, e.g., Mao & Schnei- 
der 1998, Mctcalf & Zhao 2002, Dalai & Kochanck 2002). 
With only these constraints, wc arc restricted to simple 
lens models. Circular models arc too simple because the 
lens center does not lie along the line joining the quasars. 
We therefore ask: what is the required ellipticity and posi- 
tion angle of the mass distribution? From the qualitative 
theory presented by Saha & Williams (2003), we expect 
the answer to be fairly model- independent. 

We used the "power-law" family of mass models, 

where k is the projected surface density (in units of Scrit, 
the critical density for strong lensing), ^ is the elliptical 
coordinate in the image plane, and b is the Einstein ra- 
dius. Implicit in ^ are 4 parameters: the coordinates of 
the lens center, the projected axis ratio q^, and the po- 
sition angle 9m- The exponent a determines the rotation 
curve: v{r) oc r"~^. The rotation curve is flat for a = 1, 
the isothermal case, and the physically plausible range is 
< q; < 2. With two additional parameters for the source 
coordinates, the number of parameters exceeds the number 
of constraints by one, leading to a one-dimensional fam- 
ily of models that fit the data exactly. We used software 
written by Kecton (2001) to optimize the model param- 
eters for a given choice of qm and 9m, by minimizing 

It should be noted, however, that unmodeled effects — such as a sharp 
the central third image by a large factor without significantly affecting 
2003). 



in the source plane (Kayser et al. 1989, as modified by 
Kochanek 1991). 

Figure 5 shows the contours of ™ the {qm,9m) 
plane. The line of best-fit solutions is nearly horizontal. 
The direction of the mass quadrupole is therefore well- 
constrained by the image configuration even though the 
radial density profile is not constrained. As one moves 
from left to right along this line, the mass distribution be- 
comes shallower, with a increasing from 0^2. Along the 
= curve, points arc plotted at intervals of Aa = 0.1 
with a symbol size proportional to a. The value of a is 
indicated explicitly for four representative models. The 
caustics and critical curves of those four models are shown 
in Figure 6. Over almost the entire range of a, the position 
angle of the mass model agrees with the position angle of 
the luminous galaxy disk. 

The models for which the mass and light are not quite 
co-aligned seem less plausible, on other grounds, than the 
co-aligned models. They are very shallow (a > 1.49), cor- 
responding to a rapidly rising rotation curve at a radius 
comparable to the disk scale length. They have a large 
radial critical curve (which does not exist for a < 1) and, 
within this curve, a third quasar image. These models are 
illustrated by the a = 1.7 case in Figure 6. In that case, 
the predicted third image has 7% of the flux of each bright 
image. Current radio maps limit the flux of a third im- 
age to <2% of the flux of each quasar image (Winn et al. 
2001).^" The shallow models also have small caustics (and 
consequently small lensing cross-sections), large Einstein 
radii, and large magnifications {b = 3'.'6 and jUtotai = 50, 
for a = 1.7). 

Therefore, for the most plausible range of mass distri- 
butions, the mass and light are co-aligned. Another way 
to state the result is that we can put only a very weak 
one-sided limit on the rotation curve {a < 1.49) by requir- 
ing the mass quadrupole to be aligned with the luminous 
disk (9i = 9m). If we require that the light and mass agree 
both in position angle and projected axis ratio {9i = 9m 
and qi = qm), the result is a = 0.80 ± 0.08. However, the 
latter requirement is too simplistic because it ignores the 
contribution of the dark halo, which one might expect to 
be rounder than the disk. 

Consider, for example, the case of a flat rotation curve: 
a = 1. Figure 6 shows that such a model must be rounder 
than the observed disk {qm — 0.55 > 0.41 = qi). As 
pointed out by Kecton & Kochanek (1998), the a = 1 
model can be interpreted as the projection of either an 
intrinsically flat Mestel disk, or a singular isothermal el- 
lipsoid (SIE). The disagreement of qm and qi rules out a 
pure Mestel disk, but a SIE halo is allowed. In terms of 
the three-dimensional axis ratio qsm, the projected axis 
ratio qm of an oblate ellipsoid is 

Qm = V^^lmSin^^ + cos^i, (3) 
giving q^m — 0.38 ± 0.08 in this case. Thus, assuming a 
flat rotation curve, a pure-halo model must be very oblate. 

4.2. Models with bulge, disk, and halo 

Next, we consider models with separate mass compo- 
nents for the bulge, disk, and halo. Based on the HST 

inner density cusp or central supermassive black hole — might demagnify 
the other two images (see, e.g., Mao, Witt, & Koopmans 2001; Keeton 
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Fig. 5. — Allowed power-law models of J2004-1349. Gray levels represent x < !> 4, and 9. Along the X — curve, points are plotted at 
Intervals of Aa = 0.1 with a symbol size proportional to a. Dotted lines indicate the values determined from the HST 7-band image. 



surface photometry of § 2, we described the bulge and 
disk by exponential functions. 



k(0 = i^oexp{-^/Ri), 



(4) 



fixing the scale lengths of the bulge and disk to be Rf, = 
0'.'064 and Rd = C'SS. We required the axis ratio and 
position angle of each component to be 0.41 ± 0.05 and 
— 6?5 ± 3?5, to agree with the HST-measurcd values, by 
adding appropriate terms to the x^-function with the soft- 
ware by Keeton (2001). We allowed the the central surface 
density (or, equivalently, the total mass M) of each com- 
ponent to vary. With only these two components, there 



is one degree of freedom. Even without any dark-matter 
halo, the model provides an excellent fit to the data, with 
= 0.1 and Mfc/Md = 0.16 ± 0.03 for Ax^ < 1. 
This result and its error bar are internal to our particular 
bulge/disk decomposition, but the result is fairly robust, 
as long as the bulge is taken to be much more centrally 
concentrated than the disk. For R^ = O'.'OS, the maximum 
value consistent with the HST image, the result changes 
only slightly, to M^/Md — 0.18. The other extreme, taking 
the bulge to be a point mass, gives M^/Md — 0.14. The 
latter case also demonstrates that the shape of the bulge 
is not significant; a perfectly round bulge with R^ = 0'.'064 
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Fig. 6. — Lens models of J2004-1349 for a = 1 (black), a = 0.8 (dark gray), a = 0.5 (light gray), and a = 1.7 (lightest gray). The left 
panel shows the caustics and source positions. The right panel shows the critical curves and image positions. The innermost ellipse and dot 
are the radial critical curve and third image of the a = 1.7 model. Source and image fluxes are proportional to the area of the dots. 



gives Mb/Md = 0.15. The disk scale length Rd is harder 
to measure in the HST image, because of the low sur- 
face brightness of the disk and the patchiness of the arms, 
but even for Rd — 1" the bulge-to-disk mass ratio is 0.11. 
Hence, insofar as both components can be approximated 
by exponential mass distributions, the bulge-to-disk mass 
ratio is 0.16 ±0.05. 

We used the single degree of freedom to investigate how 
this result changes with the inclusion of the dark halo. We 
described the dark halo as a spherical isothermal profile 
with a constant-density core, 

adding two extra parameters: the Einstein radius b and the 
core radius s. We do not consider flattened haloes, because 
with the present number of constraints, we would only be 
able to trace out the degeneracies between model parame- 
ters, which has already been done in detail by Koopmans, 
de Bruyn, & Jackson (1998) and Mailer et al. 2000 for 
the case of B1600-I-434. Even with a spherical halo, the 
number of parameters exceeds the number of constraints 
by one, and there is a one-dimensional family of solutions 
with = 0. We found all the solutions with 0" < 6 < 2" 
and, in each case, determined the net rotation curve and 
bulge-to-disk mass ratio. 

The spherical halo does not significantly change the de- 
rived bulge-to-disk mass ratio. The full range of values 
of Mb/Md is 0.13-0.17. The explanation is that when b 
is large enough to contribute significantly to the deflec- 
tion, the optimized core radius is also large, giving the 
halo a nearly constant density in the vicinity of the quasar 
images. This nearly constant-density sheet afi'ects only 
the optimized position and flux of the source. The rela- 



tive image positions and fluxes are controlled by the mass 
components representing the bulge and disk. 
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Fig. 7. — Rotation curve of the bulge+disk-|-halo model described 
in § 4.2. 

Although the rotation curve of the lens galaxy has not 
been measured, most disk galaxies have fairly flat rotation 
curves. From all the disk-|-bulge-|-halo models, we deter- 
mined the one with the flattest rotation curve (smallest 
mean-squared deviation from a straight line) in the range 
Rb < r < ARd. This model has Mb/Md = 0.16 and the 
rotation curve is plotted in Figure 7. The individual ro- 
tation curves of the bulge, disk, and halo are also plot- 
ted. Because the redshifts of the lens and source are un- 
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known, the velocity must be plotted in dimensionless units 
v{R)/\/2GDLT,crit, where Dl is the angular-diameter dis- 
tance to the lens, and Ecrit is the critical density for strong 
lensing. 

5. SUMMARY AND DISCUSSION 

We have shown that PMN J2004-1349 has a spiral lens 
galaxy, making it one of only a handful of spiral lenses 
known to date. Uniquely, the lens has an accurately- 
measured position, does not have a very massive neighbor, 
and is not collinear with the quasar images, which has al- 
lowed us to test whether the mass quadrupole is aligned 
with the luminous disk. The direction of the quadrupole 
is well constrained even though the radial density distri- 
bution is poorly constrained. We found that the mass and 
light are aligned within a few degrees, except for very shal- 
low mass distributions that appear physically implausible. 
This had boon shown previously for elliptical galaxies, es- 
pecially by Keeton, Kochanek, & Seljak (1997), Keeton, 
Kochanek, & Falco (1998), and Kochanek (2002), with 
increasingly large samples. However, this had not been 
shown before for spiral galaxies, owing to problems with 
those few examples of spiral lenses currently known. 

This conclusion would be weakened if tidal gravitational 
forces ( "external shear" ) from neighboring masses are pro- 
ducing some of the observed non-coUincarity. One might 
regard the close alignment of the light and mass as an ar- 
gument against a large shear. The only neighbors visible 
in the HST images, G2 and X, are faint (with fluxes -^10% 
and ~1% that of G) and are not positioned along the disk 
axis where they would produce the maximum effect. 

Using the axis ratio, position angle, and scale lengths 
measured in the HST image, we tested a model consisting 
of a bulge and disk with constant mass-to-light ratios. The 
model successfully reproduces the image configuration for 
a bulgc-to-disk mass ratio of 0.16 ± 0.05, a conclusion that 
does not change if a spherical dark-matter halo is added 
to produce a flat rotation curve. In 7-band, the bulge- 
to-disk flux ratio was found to be 0.33 ± 0.07, implying 
{M/L)b/{M/L)d = 0.5 ± 0.2 in /-band. This is a counter- 
intuitive result, because one expects disks to contain more 
young and massive stars (with smaller mass-to-light ra- 
tios) than bulges. A flattened halo would reduce the disk 
mass, but current data do not provide enough constraints 



to do more than trace out this degeneracy. It would be in- 
teresting to assess the stellar populations of the bulge and 
disk, and the possible effect of internal extinction in the 
disk, using multi-color HST images. The current y-band 
images are too shallow to provide useful color information 
for the bulge and disk separately. 

A high priority for future work is spectroscopy of the 
lens galaxy and source quasar. Knowledge of the redshifts 
of lens and source are essential for computing and inter- 
preting the mass-to-light ratios of the lens galaxy. Mea- 
surement of the circular velocity of the lens galaxy would 
break modeling degeneracies between different mass distri- 
butions with the same projected surface density (for anal- 
ogous work on elliptical galaxies, see e.g. Falco et al. 1997, 
Tonry 1998, Romanowsky & Kochanek 1999, Koopmans 
& Treu 2002). 

It may also be possible to measure the time delay be- 
tween flux variations of the quasar images, which depends 
sensitively on the radial mass profile. Or, conversely, it 
may be possible to measure the Hubble constant using the 
Refsdal (1964) method, if the radial density profile is deter- 
mined by other means. Sensitive radio observations may 
detect a third quasar image, or lensed radio jets emerging 
from the quasar cores, which can discriminate between dif- 
ferent radial mass profiles (Rusin et al. 2002; Winn, Rusin, 
& Kochanek 2003). Finally, if the sample of spiral lenses 
can be greatly expanded, it may be possible to use the en- 
semble of lens data to place statistical constraints on spiral 
galaxy structure, as has been done for elliptical galaxies 
(Rusin, Kochanek, & Keeton 2003). The enticing poten- 
tial of gravitational lensing to probe the mass distribution 
of spiral galaxies will depend on the success of at least 
some of these proposed observations. 

We are very grateful to Chuck Keeton, for his lens mod- 
eling code; Brian McLcod, for his PSF-fitting code; and 
Chris Kochanek and David Rusin, for helpful discussions 
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was provided by NASA through a grant from the Space 
Telescope Science Institute. Wc acknowledge the financial 
support of NSF grant AST-0104347 (J.N.W.), NSF grant 
AST-0206010 (RL.S.) and the Chilean grant FONDE- 
CYT/1010981 (P.B.H.). 
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